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(54) Detecting the instant of arrival of 
a seismic wave 

(57) In a seismic measuring method, a 
plurality of shocks are transmitted from 
a source (30) at the surface (32) of the 
ground to cause seismic waves to prop> 
agate through the ground to be picked- 
up by a geophone (20) located down a 
borehole (11), Traces of the signals 
received by the geophone are recorded, 
A magnitude which is a function of the 
recorded traces is calculated for each 



instant following the initial shock, and 
the statistical distribution of said mag- 
nitude at each of said instants is ex- 
amined. A threshold probability is de- 
termined, above which said distribution 
is considered to correspond to the 
existence of coherence between the 
different recorded traces for said in- 
stant. The recorded traces are ex- 
amined together in order of increasing 
time. When, for example, three con- 

(57) continued overleaf... 
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secutive instants all give a distribution 
above said threshold, it is assumed that 
the recorded signal corresponds to the 
first arrival of the shock wave at the 
geophone, and the instant of first arrival 
Is consequently determined. 
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SPECIFICATION 

Method and apparatus for detecting the instant of arrival of a seismic wave 

5 The invention relates to geophysical exploration of the subsoil, and in particular to seismic prospecting In 5 
which a shock is produced at a transmission point, thereby causing seismic waves to propagate through 
subsurface rock formations. The starting Instant of time is recorded, and the propagating waves are picked 
up at receiver points by devices such as geophones. Information about the structure of the subsoil can then 
be obtained from the variations in the received signals as a function of time, measured from the instant at 

10 which the shock producing the signals was transmitted. 

There are several methods of seismic prospecting. Some make use of waves propagating between a 
transmission point on the surface of the Earth's crust, and a plurality of receiver points which are likewise on 
the surface. Others use a transmission point that is on the surface of the Earth together with one or more 
receiving geophones located at depth in boreholes, or vice versa. 

15 In all such methods, the signals picked up by each individual geophone correspond to waves which have 15 
propagated directly from the transmission point to the receiver point, and also to waves which have followed 
a path including one or more reflections at the interface between different strata or other discontinuities in 
the make up of the subsoil. 
Since the path followed by the waves that are transmitted directly is shorter than the path of reflected 

20 waves, direct waves are always received before reflected waves. It is important to establish accurately the 20 
Instant of first arrival of a seismic wave at a receiver point in response to a given shock produced at a 
transmission point. This data is essential for determining, by suitable processing, which portions of the 
received wave trains correspond to reflected paths, and hence contain information about the stratigraphic 
structure of the region of subsoil explored. A knowledge of the instants of first arrival of the first arrival of the 

25 seismic waves is also useful for determining the transit times of sound waves in the different media through 25 
which a hole has been drilled. This knowledge can be used to calibrate the transit times measured directly in 
a borehole using an acoustic sonde. These transit times are integrated and then correlated with the transit 
times obtained from the seismic measurements in order to reset the origin of the sonic measurements and to 
rid them of the influence of borehole conditions, thereby obtaining accurate vertical profiling of the sound 

30 wave transmission properties of the formations through which the hole is drilled. 30 
A more complete understanding of these vertical seismic profiling techniques can be obtained from 
Special Publication No. 12 of the journal Geophysics which appeared in 1973 under the title "Vertical Seismic 
Profiling" by E.I. Gal'Perin. 

If seismic measurements were performed in an ideal medium that was noise-free, the signal received by a 

35 geophone after a shock had been produced at a transmission point would be flat until the first arrival of the 35 
seismic wave corresponding to the shock. At this instant, the received signal would suddenly begin a series 
of repeated oscillations. The arrival time of the corresponding seismic wave would then be the instant at 
which there was a sudden change in the slope of the signal picked-up by the geophone as it departed from a 
value of zero. On a recording of the signal from the geophone in the form of traces, called "seismic traces", 

40 the first arrival point of a seismic wave would then be easy to locate, not only by the human eye, but also by 40 
some suitable form of signal processing, since it is essential for the measuring process to be automated. One 
simple method of doing so would be to determine the average slope of the first portion of the received 
oscillation and to determine the point of intersection of a straight line of corresponding slope with the line 
corresponding to a zero value signal. In practice such a technique cannot be used because for various 

45 reasons, signals received from geophones are always affected by noise. The zero signal line is not straight, 45 
and oscillations corresponding to the reception of a seismic wave are often subject to interference. The 
methods by which the received traces are processed to detect the existence of a seismic wave arrival, and to 
determine the instant of arrival, are thus necessarily considerably more complicated. 
Further, it is common practice in seismic measurements performed between a point of the surface of the 

50 ground and a point within a borehole, to repeat the propagation measurements over successive shocks. In 50 
order to be able to stack the traces from a plurality of shock transmissions. The stacking consists of 
computing average values of the amplitudes of the traces obtained (by the simple processes of point- 
by-point summing) for each instant following a shock reference instant. This helps to improve the 
signal-to-noise ratio somewhat, not only for detecting the instant of first arrival, but also for detecting echos 

55 from transitions between adjacent strata. In present practice, when some of the traces appear to be 55 
excessively affected by noise, they are eliminated by visual inspection before stacking. Thus, these traces are 
not used in the determination of the first arrival instant of the seismic wave. Further, their information 
content is completely lost for the processing that locates strata. This means that a relatively high number of 
individual measurements need to be made to measure the conditions between a transmission point and a 

60 receiver point. Each individual measurement requires its own shock to be triggered, and a corresponding 60 
trace to be recorded. These operations take time, and the time during which the borehole installations 
remain idle adds to total prospecting costs. 

One aspect of the present invention aims to determine the arrival instant of a seismic wave after the shock 
that produced it in a manner that is more accurate and more automatic. 

65 This aspect of the invention provides a method of detecting the instant of arrival of a seismic wave 65 
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propagating through a subsoil medium between a transmission point and a receiver point, in which 
successive shocks are transmitted from the transmission point and a signal train is produced representative 
of the waves received at the receiver point during each of a plurality of observation windows corresponding 
to respective ones of said shots, a succession of identical instants in each of said observation windows being 
5 taken into consideration for each of the different signal trains and the value of the signal at each instant in 5 
said succession being recorded for the respective trains, characterised in that a magnitude which is a 
function of at least one of said values is statistically analysed to determine whether values of said magnitude 
at each instant under consideration manifests any systematic feature or coherence for the set of recorded 
signal trains, said analysis being repeated for a plurality of instants in said succession of instants, and the 

IQ instant of first arrival of the waves at the receiver point being determined as a function of the results of said 10 
analyses. The operation is preferably performed for consecutive instants in order of increasing time. 

The invention takes advantage of the fact mentioned earlier, that a plurality of signal trains, each 
corresponding to a respective transmission of shock waves, is recorded from each receiver point to be 
subsequently combined in order to improve the signal-to-noise ratio, particularly when at least one of the 

15 transmission and receiver points is located in a borehole. Instead of determining the instant of first amval 15 
from the stack of said recorded trains, all the information contained in each trace is taken into account when 
determining said Instant of first arrival relative to the instant the shock was transmitted. 

In one implementation of the invention, a minimum threshold is fixed in advance for the probability of the 
distribution of said magnitude at any given instant corresponding to coherence in said train. It can then be 

20 decided that said probability threshold being exceeded for one instant, or for a combination of instants such 20 
as three in succession, is indicative of the first arrival of the received seismic wave. The instant of wave 
arrival can then be fixed as a function of the Instants in said combination. For example, the instant of first 
arrival may be chosen as the instant preceding a series of three consecutive Instants in increasing time order 
during which the chosen probability threshold is exceeded. 

25 In a particularly advantageous implementation of the invention, when it is desired to determine the instant 25 
of first arrival of the seismic wave in real time, said instant is sought on the basis of some minimum number 
of traces which are analysed using the above-defined method. If such analysis fails to determine an Instant of 
first arrival with the required probability, the newly acquired trace is added to the previously analysed traces, 
and the procedure is repeated. 

30 The magnitude whose distribution isanalysed may simply be the value of each of theseismictracesatthe 30 
instant under consideration. It is also possible to use a function of said value and at least one other value in 
the train. For example, we have been successful using a combination in which the magnitude taken into 
consideration at each instant is a function of the average slope of the train in the vicinity of the instant under 
consideration. 

35 Using the method outlined above, it is possible to take advantage from all of the seismic traces, even the 35 
noisy ones, recorded for propagation between a given transmission point and a given receiver point. There 
is no need for prior visual inspection. The method is systematic. It does not make use of techniques which 
simulate the kind of analyses which a human observer might perform: techniques that run into difficulties 
when the traces are masked by a relatively high level of noise. Finally, it can be performed quickly. 
40 The invention also provides apparatus for detecting the first arrival of a seismic wave, said apparatus 40 
comprising a combination of means which are organised to perform the method outlined above. 

Further explanation, and a description of particular implementations of the invention are given below by 
way of example with reference to the accompanying drawings, in which: 
Figure 1 is a diagram showing a borehole seismic prospecting installation; 
45 Figures 2A and 2B are diagrams of seismic traces for illustrating the way the invention is implemented; 45 
Figure 3 is a table diagrammatically representing a series of seismic traces in the form of samples taken at 
discrete time intervals; 

Figure 4 is an outline of a table for determining a probability for use in the method in accordance with one 
aspect of the invention; and 

50 Figure 5 Is a flow chart outlining the operation of a data processor when seeking a first arrival time in 50 
accordance with the invention. 

A seismic borehole probe 10 is lowered down a borehole 11 (Figure 1) on the end of a cable 1 2 which 
serves to move the probe along the borehole and to connect it electrically to a mobile laboratory 14 on the 
surface. The cable passes over a sheave 15 at the top of a drilling rig 16 and drives a pulley wheel 18 
55 tangentially in such a manner that the number of pulley wheel revolutions is proportional to the depth of the 55 
sonde 1 0. In the position shown, a geophone 20 has a sensitive face located in or close to the outside wall of 
the tool 10. This face is pressed against the inside wall 21 of the borehole 11 by means of two anchoring 
members 22 which are fixed to the tool 10 and which press against the opposite face 23 of the borehole. 
A source of seismic signals 30 Is located in a water-filled trench 33 in the surface of the ground at a suitable 
60 distance from the borehole 1 1 . The source is, in theory, shaped so as to produce a high powered shock of 60 
short duration in the ground, and may be of any well known form, eg. an air gun. 

The source 30 is controlled from the laboratory 14. A hydrophone 34 placed close to the source 30 detects 
the instant at which the wave departs into the ground. This instant serves as a reference point for each shot, 
i.e. the traces picked up by the geophone 20 are all timed from the moment the corresponding shock 
55 produced by the source began to propagate. 
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Seismic waves propagating between the source 30 and the geophone 20 can follow multiple paths. The 
direct path is referenced 35, A path with two reflections, one at a stratum 36 and another at the surface 32 of 
the ground is referenced 37. It results in a wave that is travelling downwards when it reaches the geophone 
20. A third path 38 has a single reflection at a stratum 39, and results in a wave travelling upwardly when it 
5 reachesthe geophone 20. 5 

Because the paths 35, 37 and 38 are of different lengths, the electric signals transmitted by the geophone 
to the laboratory 14 are shifted in time relative to the initial pulse. Traces 40, 41 , 42, 43, and 44 in Figure 2A 
are records of the amplitudes of the oscillations picked up by the geophone 20 during five successive shots 
from the source 30, plotted as a function of time TV measured from the initial reference instants which are off 
10 to the left of the diagram. The pauses between successive shots are long enough to ensure that each io 
recorded trace corresponds to substantially all seismic propagation paths capable of supplying measurable 
data. 

The left hand point portion 50 of the trace 44 is practically flat and parallel with the time axis up to a point 
51 corresponding to an instant TT% where the signal Is deflected downwards (trough 52) before climbing to a 

15 peak 53 at a greater amplitude than the level 50. The set of curves 40 to 44 are manifestly similar to one 15 
anotherto the right of the instant TTi. The similarities between the curves are representative of the existence 
of systematic features or coherence in the propagation of seismic waves between the transmission point and 
the receiver point regardless of the shot that produced the waves. This coherence between the traces is 
manifest, inter alia, from the way the points of first arrival at the instants TTi can be seen by eye to be in 

20 alignment. However, it should be observed in the case of curve 40 whose left portion is degraded by noise, 20 
that it is not easy to locate the exact instant of first arrival by eye or by any method that is functionally similar 
to the human eye and brain. 

The signals picked up for each trace (via Input 63) are acquired and digitised under the control of a control 
and processing unit 64 in the laboratory 14. The control and processing unit 64 receives depth data in the 

25 form of pulses from the pulley 18 via a link 65. Unit 64 is connected to the source 30 and to the hydrophone 25 
34 via respective conductors 62 and 66. A magnetic recording unit 67 stores digitized data corresponding to 
each reco rded trace. 

Let there be n traces Ti to Tn acquired under the control of the device 64 from a sequence of n shots (Figure 
3). The amplitude of each trace is represented by digital values Xi(k) obtained by sampling over regular time 

30 intervals at instants A:, where /r lies in the range 1 to /. For each trace the instants 1 to fare determined from an 30 
origin constituted by the instant at which the corresponding shot, was generated. Values having the same 
index A:/ such as Xi{k), X2(k), . . . Xn(k), correspond, to the same transit times for the seismic wave between the 
source 30 and the receiver point or geophone 20. Nonetheless, these values are not identical because of the 
noise to which the signals are subject. 

35 Signals picked up at instants /c prior to the instant of first arrival of waves at the geophone 20 can only 35 
correspond to noise Nj(k) for each trace /. On the basis of this observation, each group of values concerning a 
particular instant Ic, {for Ic in the range 1 to fl, can itself be considered as constituting a sample Hk taken from a 
set of numbers that obey different laws of statistical distribution depending on whether the sampling instant 
(c is before or after the instant of first arrival. 

40 Supposing that the noise is Gaussian, the mathematical expectation associated with each sample Hk prior 40 
to the first arrival instant must be zero. In contrast, the expectation must be non-zero for the instants which 
Immediately follow the first arrival. 

Since the number of recorded traces may be relatively small, the distribution of each sample Hk 
corresponding to an instant Ar is supposed to obey Student and Fisher's law. 

45 Student's t-variable is calculated for each sample of order Arusing the equation: 45 

n 

2X|(k)2(n-1>2 

50 I 50 



(1) V = 



n n 
n2Xi(k)2-(2X|(k))^ 
I I 



55 



55 



in which X;(k) is a function of the value Xi(k) and possibly other values, as is explained further on. 
According to Fischer's law, the variable tk has a probability distribution given by the equation : 

60 1 r((no+l)/2) * 

(2) p(tk) = 

^/2?k r {nJ2){^ +t2/no)((no+1 )/2) 
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in which the function r obeys the relationship Hn+I) = nHn) and Pq is the number of degrees of freedom. 

When there are n independent seismic traces, the number of degrees of freedom, or number of draws or 
samples that can be drawn from a set of Gaussian distribution, is (n-1). 

As n increases, the distribution of the variable tk gets closer to a normal or Gaussian distribution 
5 corresponding to the distribution of the set from which the draws are taken. When the number of items in 5 
each sample under consideration (ie. the number of traces) exceeds 30, the t distribution approximates very 
closely the normal distribution. 

For smaller numbers of traces, the value of the probability p(tk) may betaken from a Student and Fischer 
table, instead of being calculated from equation (2). It is then possible to obtain directly or by extrapolating 
10 the probability Pk of the set from which the sample was taken having an average value of zero {only noise ig 
present) and the probability (1 -pk) of the average value being other than zero, ie. of the set from which the 
sample was taken corresponding to a superposition of noise and signal which has a similar or coherent 
effect in each of the traces. 

Figure 4 outlines a table 83 of values of the variable f which may be stored in a memory of the apparatus 68 
15 for detecting the instant of first arrival. The table has two sets of indices for defining rows and columns 15 
respectively. Each row corresponds to a particular number of degrees of freedom Rq. Each coium 
corresponds to a given probability p that a sample is taken from a set having an average value of zero, or to 
the probability (1-p) that the reciprocal condition istrue. 
Thus, given a set of traces T^ to Tn, which have been amplitude-sampled at discrete instants Ar in the range 
20 1 to /(Figure 3), Student's t-variable t^ may be calculated for each instant Ar. The table of Figure 4 provides, for 20 
the number of degrees of freedom n-1 and the value tk in that row, the corresponding probability (1 -pk) that 
the sample Hk constituting values Xi(k) to Xn(k) was taken from a set having a non-zero average value and 
thus indicates the presence of a useful signal. 
Consequently, the samples Hk may be examined in order of increasing time until the probability (1 -pk) 
25 exceeds a predetermined threshold, eg. 0.99. To make the operation more certain, it may be decided, as in 25 
this example, that the presence of a coherence signifies the existence of a useful signal only after three 
consecutive samples Hk In time increasing order have all given a probability (1 -pk) which exceeds said 
threshold. 

The instant of first arrival is then defined as the instant immediately preceding the three consecutive 
30 samples which satisfied the required condition. 30 

Curve 80 in Figure 2A shows the value of the probability (1 -p) obtained using the above method on the 
five traces 40 to 44. By reading the curve 80 from left to right, it can be observed that the instant of first arrival 
TTi is the first instant after which the probability (1 -p) is stable at its maximum value for several consecutive 
instants of time TT. 

35 In practice, the sample Hk may comprise variables which, although dependent on the instant k under 35 
consideration in the trace detection interval, are not necessarily the values or the amplitudes of the traces at 
said instant /r. 

Thus, in general terms, the sample Hk comprises the values of a variable Xi(k) derived from each trace Ti. 
For example, instead of using the amplitude X||k) of the trace T| as the value for Xi(k), the difference between 
40 two consecutively sampled values of the trace could be used, ie. Xi(k) - Xi{k-1). It has also been observed 40 
that good results can be obtained by letting X-,(k) have the value of the gradient of the least squares fit line 
through five consecutive samples from the instant k, ie. through samples X||k) to equation has been used: 



45 (3) Xi(k) =^ [- 2xi(k) - Xi(k+1 ) + Xi(k+3) + 2X||k+4)l 45 

The variable derived from equation (3) has the advantage of eliminating any DC component that may be 
present in the seismic traces. 

50 Figure 2B shows traces 140, 141, 142, 143 and 144 which have the same information content concerning 50 
seismic propagation between the points 30 and 20 as the traces 40 to 44 in Figure 2A, but each of the traces in 
Figure 28 has had additional Gaussian noise added thereto. 

A curve 81 shows the resulting values of the probability (1 -p) obtained using the method in accordance 
with the invention. It should be observed that the curve 81 still enables the point TTi to be accurately 
55 determined, even though determining its position from any one of the traces 140 to 144 has been made 55 
difficult, if not impossible. 

Figure 5 summarises the operation of the control unit 64 shown in Figurel for detecting the first arrival 
time as the seismic traces are being acquired for waves propagating between a transmission point 30 and a 
receiver point 20. It is assumed that the search for the instant of first arrival is commenced once some 
60 n^inimum number of traces, eg. 5, has been recorded. 60 
After initialisation (block 200), a counter for counting the time instants k is set to unity and the value of the 
variable X is calculated for the trace n currently being acquired, ie. Xn(k), see block 204. The value 
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n~1 

t 

5 is taken from a memory 206 which is used for storing the sum of the values of the variable X at instant Ar as 5 
calculated for all the traces obtained so far. The stored sum is added to the newly calculated value from block 
204 to obtain (block 208) a new sum value 

10 n 10 

2Xi(k). 
I 

The value Xn^(k) Is then calculated In block 210. Then the sum of said squares 
15 15 

n-1 

2Xi2(k) 

I 

20 20 
as calculated so far is taken from a register 21 2 and a new sum 

n 

25 SXi^fk) 25 

1 

is calculated in a block 214 from the new value provided by block 210. 
A check is then made to see whether the number of traces acquired so far is greater than the minimum N 
30 required for investigating the distribution (block 216). If not, eg. \fn is less than 5, a check is made to see if all 30 
the instants k of an observation window have been dealt with (block 220), and the system returns to block 
204 to repeat the procedure outlined above. This is continued with all /instants in the interval have been 
dealt with. 

When this happens, (k=f in block 21 8), the number of the trace is incremented (n=n+1, see block 222) and 
35 the next trace is acquired by returning to block 202. 35 
When the number of traces already acquired is greater than or equal to the minimum N (answer YES in 
block216). Student's variable (t^) is calculated in block 224 in the manner already described. 

The corresponding probability (1-Pk) is then looked up in the Student-Fisher table (block 226) described 
with reference to Figure 4. Naturally, the value of (1 -Pk) may be interpolated should the calculated value of t^ 
40 turn out to fall between two adjacent stored values in the row corresponding to n-1 degrees of freedom. 40 
A check is then made to see whether k is greater than or equal to 4 (block 228), and if not, k is incremented 
by unity and the system returns to block 204. 

Otherwise, a test is performed to see whether the probabilities obtained for instants k, k~ 1 , and k-2 are all 
above a threshold probability Ps which is predetermined and stored in an external register 230, eg. having a 
45 value of 99%. If this condition is satisfied, (answer YES to block 232), that means that the signal 45 
corresponding to the first seismic wave arrival has been found by reading through all n traces together in 
order of increasing time. The instant of first arrival is then taken to be k~3 as indicated in block 234. 

In the event that the condition tested in block 232 gives the answer "NO", a test is made to see whether the 
entire measurement window for the trace n has been examined (block 236). If not, k is incremented by unity 
50 (block 238) and the system returns to block 204. Otherwise, if all the samples from trace n have been 50 
examined without it being possible to determine the instant of first arrival, the system moves on to the 
acquisition or reading of the next trace and returns via block 222 to block 202, and the process is repeated. 

Once enough traces n have been obtained for finding the first arrival time with sufficient accuracy, and the 
arrival time has actually been determined (block 234), the process is terminated. 
55 The above operations can be implemented with the automatic processing equipment usually to be found 55 
in mobile logging units. The only specific requirements are two buffer registers (blocks 206 and 212) for the 
sums and a three position memory for the values (1 -pj* 

Instead of doing such processing in real time, the processing can be done in deferred time, provided a 
plurality of previously recorded traces is available and stpred in the magnetic memory 67. The number 
60 Nm nriust be much greater than the minimum number N. The operations outlined in the set of blocks 204 to 60 
226 is then performed for each instant in increasing time order until three values are found which satisfy 
the condition in block 232. 
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CLAIMS 

1 A method of detecting the instant of arrival of a seismic wave propagating through a subsoil medium 
between a transmission point and a receiver point, in particular where at least one of said points is located 

s underground and in which successive shots are transmitted from the transmission point and a signal tram 5 
(i) is produced representative of the waves received at the receiver point during each of a plurality of 
observation windows corresponding to respective ones of said shots, a succession of instants |k) counting 
from a reference instant in each of said observation windows being taken into consideration said succession 
being identical for all the trains, and with the value Xilk) of each signal train at each instant (k) being recorded, 

in forthe respective train, wherein a magnitude Xi(k) which is a function of at least one of said values Xi(k) is 10 
statistically analysed to determine whether the distribution of the samples (H^) constituted by the values of 
said magnitude which correspond to an Instant k manifests any systematic feature or coherence of the set of 
trains received, said analysis being repeated for a plurality of successive instants k. and the instant of first 
arrival of the waves at the receiver point being determined as a function of the results of said analysis. 

2 A method according to claim 1, characterised in that the statistical analysis operation IS repeated for 15 

said successive instants in the order of increasing time and said instant of first arrival is determined as a 



function of at least successive instants. u i^ r*. *i. --k^w \ i» avoh 

3 A method according to claim 1 or 2, characterised in that a minimum probability threshold (Ps) is fixed 
above which said distribution of the sample (Hk) at an instant k is considered to correspond to the pressure of 



20 



20 coherenoB^^^^^ according to claim 3. characterised in that the received signals are considered to correspond 
to the first arrival of the seismic wave once said coherence has been detected for a predetermined 
combination of a plurality of successive instants (k-2,k-1,k). 
5. A method according to 4, characterised in that the instant of first arrival is taken to be at a 
25 predetermined instant (k-3) preceding said combination of instants. ^ • , x^^^ ^ 

6 A method according to any preceding claim, characterised in that the instant of first arrival of the 
seismic wave is sought on the basis of a minimum number of trains which are analysed at instants in order 
of increasing time in the observation window, and, in the event of an unsuccessful search, in that the search 
isrepeatedbyaddingatleastonenewlyacquiredtraintothesetoftrainsunderconsideration. 
30 7. A method according to any preceding claim, characterised in thatthe presence of coherence In the 30 
signal trains (Ti to T„) is determined by means of the Student-Fisher test. . ^ . 

8 Amethodaccordingtoanyprecedingclaim,characterisedinthatthemagnitudeXi{Wtakenmto 

consideration at each instant (k) is a function of the valuexi{k) of the signal in each train at said instant and of 
at least one othervalueofthe signal in said train. . ^ ^,,y.^t. 

« 9 A method according to any preceding claim, characterised in that the magnitude Xi(k) taken into 35 

consideration at each instant (k) is a function of a feature associated with the slope or the average slope of 
the signal in each train (i) in the neighbortiood of the instant under consideration (k). 

10 Apparatus for detecting the instant of arrival of a seismic wave propagating through a subsoil 
medium from a transmission point to a receiver point, the apparatus being of the type th^;i'"'='"^f ^-T^^^^^^^ 
40 for acquiring picked-up during an observation window following the transmission of a shock, and means for 
^ recordingsaidacquiredsignalsinthefomioftracesofapluralityofsuccessiveshocks characten^^^^^^^ 
i,e apparatus further includes means for examining the distribution (H.) of a magnitude ojtamed from the 
set of recorded traces at each instant in the observation window, means for determining the P^f ^^^i^ /Jj^ 
the distribution of said set of traces manifests the existence of coherence in the set "^''.^^ Jf/^^^"^^ 
45 corresponding instant, and meansfor detecting when said probability excewis a predetermined threshold 
1ApparLsaccordingtoclaim10,characterisedinthatitfurtherincludesareg.ster^^ 
values of said magnitude, and a register for summing the squares of the values of said magnitude, and a 
memory for storing a probability distribution table for a Student-Fisher variable. 
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